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UNCLASSIFIED 


Nonlinear  response  of  a superlattice 

Fla,  T.,  Birnir,  B.,  and  Allen  Jr.,  S.J. 


Abstract — In  this  work  we  study  the  nonlinear 
dc  and  ac  response  of  a superlattice.  We  de- 
velop asymptotic  theories  for  the  I-V  curve  and 
ac  gain  valid  in  both  the  very  low  frequency  qua- 
sistatic limit,  the  intermediate  frequency  limit 
and  high  frequency  limit.  We  also  explore  the  I- 
V curve  and  ac  gain  for  a quite  extensive  param- 
eter regime.  Parameter  regimes  which  have  gain 
and  positive  differential  conductivity  are  identi- 
fied for  sufficiently  high  ac  amplitudes.  These 
parameter  regimes  will  be  important  to  under- 
stand the  nonlinear  state  of  solid  state  sources  in 
the  THz  range. 


I.  Introduction 

We  will  in  this  work  report  on  nonlinear  ac  and  dc  con- 
ductivities in  a superlattice  favorable  for  a THz  source 
[1].  Esaki  and  Tsu  [1]  proposed  to  use  a semiconduc- 
tor superlattice  to  to  achieve  a device  with  a current- 
voltage  characteristic  with  negative  differential  conduc- 
tivity (NDC)  above  a certain  critical  dc-bias  propor- 
tional to  the  geometric  mean  collision  frequency  (z/p). 
Subsequently,  Ktiterov  et  ah  [2]  gave  a simple  model  for 
the  miniband  electron  distribution  function  for  a super- 
lattice including  both  momentum  and  energy  relaxation. 
They  showed  that  for  dc-bias  in  the  NDC-regime  and 
linear  ac-drive  that  there  was  ac-gain  for  frequencies 
slightly  below  the  resonant  Bloch  frequency.  We  will 
further  investigate  Ktiterov’s  theory  into  the  nonlinear 
ac-drive  regime.  It  turns  out  that  the  problem  of  calcu- 
lating currents  in  Ktiterov’s  Boltzmann  model,  can  be 
reduced  to  studying  the  momentum  equation  and  en- 


ergy equation [3,  4]  .The  electron  density  equation  for 
homogenous  E-field  is  simply  p0  = const  . 

The  model  we  will  concentrate  on  is  based  on  the  fluid 
moments  of  the  Ktiterov’s  Boltzmann  equation,  but  in 
a nondimensional  form  where  the  essential  parameters 
of  the  theory  is  stressed  from  the  beginning.  This  em- 
phasis has  simplifying  implications  for  scaling  of  theo- 
retical, simulation  and  experimental  results  and  param- 
eters. Our  goal  is  to  find  conditions  under  which  a dc- 
biased  and  nonlinearly  ac-driven  superlattice  can  have 
both  gain  and  positive  differential  conductivity.  The  co- 
incidence of  the  regime  of  gain  of  small  amplitude  oscil- 
lations and  negative  differential  conductivity  imply  that 
domain  formation  is  an  obstacle  to  the  construction  of  a 
superlattice  THz  source.  The  ac  E-field  is  assumed  to 
give  rise  to  such  high  electric  fields  in  the  medium  that 
the  dc  and  ac  conductivities  in  the  superlattice  is  non- 
linearly modified  . We  want  to  explore  conditions  under 
which  the  dc  conductivity  is  locally  of  the  positive  dif- 
ferential type  and  at  the  same  time  that  the  nonlinear 
conductivity  give  rise  to  gain.  The  model  we  choose  to 
investigate  is  the  momentum  p^v  and  energy  p0£  fluid 
moments  of  a Boltzmann  equation  with  a simple  energy 
and  momentum  collision  operator  . We  find  the  follow- 
ing equations  in  agreement  with  several  earlier  works  : 


fd  , ed2E , 
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By  introducing  new  variables  v€  = (e  — €to )d/K  , 

v®  = — £To)d/h  = (Aid/2h)  T][2kS!{)\  and  defin- 

es kBT0) 

ing  the  Bloch  frequency  ojb  = edE/H  , we  find  the 
alternative  description 


{-Qt  + vp)vp  + WBVe  = UBV°e,  (2) 

7^  + Ve)ve  - UBV  - 0 (3) 

These  equations  can  now  be  put  into  a nondimensional 
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form  by  scaling  the  velocities  by  v®  ( alternatively  by 
(Ljbo  or  duo) and  time  and  frequency  by  the  background 
Blochfrequency  uoBo  = ~ (alternatively  uo  or  the  the 
geometric  mean  collision  frequency  which  we  also  will 
discuss  below  ).  We  then  find  the  following  nondimen- 
sional  form  of  the  equations  with  frequencies  and  time 
scaled  with  respect  to  the  dc  Bloch  frequency.  Here 
we  have  by  abuse  of  notation  used  the  same  symbols 
for  the  variables  uoB  ~ 5 ► 
t -»•  tu>BO,  ( V£,V ) -»  (v£,v)/v£0 


,d 

+ l’v)v  + ubv£  = wB) 

(4) 

{-Qt  + "eH  - UBV  = 0. 

(5) 

An  alternative  useful  scaling  of  frequencies  and  times 
of  this  equation  are  given  in  terms  of  units  of  the  ge- 
ometric mean  collision  frequency  (c.f.  discussion  be- 
low) ug  = (z^z^)0,5  which  means  t — > t = tug; 

Ve,p  — > k£P,-^;ujbo  ->w0  = ™;wB-tw B ~ Tf 

The  dc-solutions  of  these  equations  for  a given  con- 
stant E- field  £o  = Vo fL  corresponding  to  a constant 
Bloch  frequency  to  bo  is  given  in  these  scaled  variables 
by 


( I'd  = (vp  - Ue)/2  = vd/u 0,  I'd  = (1/ Kep  - 
Kep)/2  ~ 0.35)  which  for  high  enough  Wo  is  close  to 
the  normalized  Bloch  frequency  1.  For  fixed  Bloch  fre- 
quency (u/o)  or  fixed  frequency  (d>)  this  resonance  trans- 
late to  Qr  = ±(o;2  - ^d)0'5  or  u0r  = ±(u)2  — f>2)0'5 
respectively.  A short  calculation  will  show  that  there  is 
linear  gain,  i.e.  <T\  = Re((J /)  < 0 only  if  coq  > 1 and 
in  the  frequency  regime  to  < toz  < tooc  = 1.  Here 
OJ^near  to  tor)  is  the  crossover  frequency  from  gain  to 
no  gain.  Consequently  there  is  no  coincidence  between 
gain  and  PDC  in  the  linear  limit  tO\  — > 0 We  remark 
that  no  resonance  is  found  for  Ud(—  -35/uJq  )>  1 , i.e. 
u)q  < I'd  or  do  < Vd  respectively.  If  we  assume  that  the 
ratio  of  the  collision  frequencies  are  reasonably  fixed,  we 
see  that  the  description  of  the  dc  and  ac  response  can  be 
given  in  terms  of  the  three  nondimensional  parameters 
CC\  = tO\tOo  , Co  = ojooq  and  loq  which  is  nothing  but 
the  ac-drive  amplitude,  frequency  and  Bloch  frequency 
measured  in  units  of  the  geometric  mean  collision  fre- 
quency. These  three  variables  are  also  reasonable  to  use 
since  for  given  geometric  mean  collision  frequency  they 
correspond  to  fixed  frequency  and  ac-amplitude  as  we 
vary  the  dc  bias. 

II.  AC  AND  DC  CONDUCTIVITIES  WITH  A 
NONLINEAR  AC-DRIVE 


V0=  KepOJo/(  1 + LOq)  = <T0W0;  Ve0=  Wo/ (1  + Wq) 

(6) 

Here  we  have  introduced  the  scaled  dc  conductivity 
(7$  . We  have  also  introduced  the  relative  Bloch  fre- 
quency in  units  of  the  geometric  mean  collision  frequency 
too  which  is  normalized  to  1 at  the  maximum  of  the  I- 
V curve.  Another  quantity  we  will  be  interested  in  is 
the  linear  response  of  the  velocity  and  energy  due  to  a 
harmonic  perturbation  in  the  E-field  drive  correspond- 
ing to  00  b ~ 1 +eOOBl(t),  C OBl(t)  = — 0/iCOs(cJt)  = 
— exp(— iuot)  + c.c..  e is  just  an  ordering  parameter. 

We  now  find  the  0(e)  scaled  ac- response  (the  dc  and 
ac  response  for  a nonlinear  ac  drive  is  too  long  to  be 
included  explicitly)  in  agreement  with  Khiterov’s  model 
ac-conductivity 


_ , Vl  S_  WjWo  <7i(w)  s 

Vl  " (0«i  )— ] 

/ &i  \ £o 

' aei  ' w§  + 1 — (ww0)2  - iu>u>o(nEp  + l//tep) 

/ 1 — W§  - iulOjQ  / nep  s 

^ oj0(2/k£P  - iww0) 


(7) 

(8) 
(9) 


w2 


v€  — ito , vp  = Up  — ito  (10) 


We  immediately  see  that  there  is  a resonance  in  the 
ac-conductivity  at  frequency  tor  = ±(1  — I/2)1/2 


We  summarize  our  theoretical  (there  is  no  space  to  give 
the  theory)  and  numerical  investigation  below(see  some 
of  the  results  of  I-V  curves  and  ac  response  in  Fig,  2a,b  ) 
This  has  eventually  to  be  substantiated  by  experimental 
results. 

1. The  solution  of  our  model  in  the  high  frequency  limit 
approaches  a limit  cycle  with  a timescale  r — 0(1 /v) 
and  with  leading  order  oscillation  frequency  and  renor- 
malized resonance  given  as  coi  = (I  — Vd )1//2  if  & d — 
Vd  Jo(2uq/w)  < 1.  For  uq  > 0(1)  there  is  a narrow 
layer  of  PDC  and  gain  in  the  domain  u;o  6 [to  — 6>Co] 
where  S = o(l)  and  similar  indications  near  the  sec- 
ond resonance  ^ 2 Co  for  £q  1.  There  is  a change 
in  signature  of  the  I-V  curve  for  small  too  and  near 
Cui/Co  = 2.405  (Jo(^)  =0).  The  corrections  of  0(vd) 
is  quantitatively  important  in  our  theory  in  contrast  to 
many  authors  who  use  a model  with  only  one  collision 
frequency. 

2.  The  quasistatic  theory  and  intermediate  frequency 
theory  to  < O(ud)  is  characterized  by  a decay  with 
timescale  T (no  oscillations)  to  the  limit  cycle,  a linear 
shift  of  the  maximum  of  the  I-V  curve  u)qc  ~ u)\  and  a 
small  region  of  coinciding  gain  and  PDC  for  U)\  > 1. 

The  limit  cycle  can  be  described  in  a threedimen- 
sional phase  space  plot  of  {u7#(£)}  v(£)}.  Since  we  will 
basically  be  interested  in  conductivities,  we  will  study 
the  projection  of  the  limit  cycle  on  the  Blochfrequency- 
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velocity  plane,  ^5  (i) , } .The  limit  cycle  will  have 

a basic  period  of  T — 2'k/u,  but  can  have  a complex 
behavior  with  one  or  several  loops  at  high  frequencies 
u>  1 and  will  due  to  the  stiffness  of  the  system  have  a 
tendency  to  align  along  the  I-V-  curve  in  the  quasistatic 
limit  for  small  Q . 


A.  I-V  curves 


T 

The  averaged  velocity  Vq  = ^ f vdt  over  a basic  period 

0 

T = 2'k/u)  of  the  limit  cycle  is  the  dc-  component  of 
the  periodic  velocity.  DC  I-V  curves  (see  examples  in 
Fig.  2a)  in  this  context  is  nothing  else  than  a graph  of 
vo(u>q\QuQ)  as  a function  of  the  Blochfrequency  with 
the  amplitude  and  frequency  of  the  drive  as  parame- 
ters.To  measure  this  property  for  a nonlinearly  driven 
superlattice  numerically  one  has  to  wait  until  all  tran- 
sients has  died  out  such  that  the  solution  sticks  to  the 
limit  cycle.  For  the  parameter  regime  and  the  collision 
frequencies  we  are  interested  in,  it  is  suffient  to  wait  4-5 
periods  before  one  do  the  averaging  around  one  period 
of  the  limit  cycle.  The  collision  frequency  for  a GaAs 
superlattice  has  been  estimated  by  a fit  to  I-V  mea- 
surements fo  a specific  sample  given  on  the  Quantum 
Institute  homepage  (J.  Scott)  to  be  \jv£  ~ 0.3ps  and 
1 jvp  ~ 0.15ps  . For  a InAs/AlSb  superlattice  a rea- 
sonable guess  for  collision  frequencies  are  \/v£  — 0 6ps 
and  1 jvp  ~ 0.2  — QAps  . Therefore  estimates  of  the 
collision  frequency  ratio  are  Kep(GaAs)  ~ 0.707  and 
Kep  £ [0.6, 0.8]  . The  geometric  mean  of  the  colli- 
sion frequencies  are  therefore  estimated  as  vg  ( GaA$ ) ~ 
4.7 THz  and  vg(InAs/AlSb)  £ [2,3 ]THz.  With 
these  parameters  given,  the  Bloch  frequency  can  be  de- 
duced, i.e,  u)Qy  u)j  u?i  = 1 corresponds  to  the  physical 
frequencies  — vg  To  help  the  reader  we 

give  the  following  table  which  scale  our  nondimensional 
theoretical  and  numerical  results  to  different  type  of  su- 
perlattice  experimental  settings. 


4.7w0(or  Wi)(  J 

V0(U) 

0-294(^)J^^(or^i)eV  1 

r i/fjUJ 

J ~ 2tt 

0-748 (r^pTHz 

VeO 

2^  ii>(«)V9  X 10e?n/S 

K — Al 
^ “ 2 knT 

0 3l68(^)(3o^FT 

(7 

( N \-l  no  S Lrnf 

vlOO'  10 u>cm  3 20 0.1615^i7n 

Tub. I Seating  of  some  physical  quantities  from  lion  dim  cusioual  ones. 

IIcrf;  a = a00.393mho  „„<i  a = " For  scftli„K 

of  physical  (7  ^p  with  nondimensional  (T  J interchange  (7  ^ (7  \ and 

^0  — 3 ► CTlp- 

Here  F(k)  = is  plotted  in  Fig.l  and 

can  also  be  used  to  infer  the  nonlinear  temperature 
dependence  of  the  dc-bias  and  ac-field  given  data  for 
the  nonlinear  mean  energy  through  our  simulation  (we 


Figure  1 Here  we  have  plotted  the  function 

Ji(k)//0(k)  where  the  nondimnnsiounl  parameter  AC  — 
scaling  of  this  pnmmctcr  for  different  supcrlatticcs  is  given  ii 


F(k)  = 

A , 

2k bT  ' T1,c 


do  not  discuss  the  energy /temperature  results  here)  or 
by  experimental  means.  v£jy(u)iu)iiu)Q)  =.^^F{k)  => 
T{w,Qu^o)  = ^F-\-^-dvED). 

This  nonlinear  temperature  might  be  an  important 
parameter  to  determine  selfconsistent  models  of  collision 
frequencies,  especially  for  temperature  regimes  where  ki- 
netic energy  is  comparable  to  the  phonon  energy.  We 
have  also  defined  a characteristic  collision  length  Lrnf  = 
Vso/Vg  for  superlattices.  The  dc  current  given  for  a 
nondimensional  conductivity  <7o  in  a superlattice  with 
electron  density  no,  transversal  crosssection  of  the  mesa 
S and  length  L — Nd  now  scales  as  a function  of  dc 
voltage  Vo,  ac  electric  field  E\  = V\/L  and  frequency 
/ by  Idc  — aVo  ([/]  — Ampere  if  [Vo]  — Volt).  The 
ac-current  and  conductivity  scale  the  same  way  given  an 
ac  nondimensional  conductivity. 

B.  Nonlinear  ac  response 

We  will  now  measure  the  nonlinear  ac  response  (see 
examples  in  Fig.  2b)  of  the  superlattice  velocity  and 
energy  at  the  first  harmonic  of  the  limit  cycle.  The 
first  harmonic  can  be  determined  directly  through  a dis- 
crete fouriertransform  of  velocity  and  energy  data  from 
a complete  period  of  the  limit  cycle.  Alternatively,  the 
real  part  or  cosine  component  of  the  first  harmonic  can 
be  determined  by  averaging  the  superlattice  equations 
over  one  period  assuming  again  that  we  are  at  the 
limit  cycle  such  that  periodicity  can  be  assumed,  i.e. 
u(£)  — v{t  + T).  This  leads  to  the  averaged  equations 


vvvD  + veD  - ~2Vic 

= 1, 

(11) 

VeVeD  ~VD  + — V£lc 

= 0. 

(12) 

The  interaction  terms  T^^iCjelc  has  an  inter- 
esting physical  interpretation  as  the  rescaled  averaged 
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Figure  2 In  Fig.  2a  we  hnve  plotted  the  normalized  dc  I-V  curve 
pnru m ctcrized  by  the  velocity  nveraged  over  one  period  of  the  limit  cycle  as 
a function  of  normalized  Bloch  frequency  £ [0,  10]  and  for  a fixed 
nonlinear  ac-drlvc  amplitude  5*0.  These  normalized  I*V  curves  can 

be  scaled  to  I-V  curves  for  different  experimental  parameters  by  the  scaling 
laws  given  in  Tab.l  and  the  relation  Vq  = (JOJQ.  The  curves  arc  labelled 
ns:  dnsh-dotted  curve  ^ LO  = 0*1,  curve  of  + ’s  ‘ ^ “ 1.0)  curve 

of  o’s  — ♦ CJ  = 2.079  •—  2.405)and  solid  curve  ^ U)  — 5.0* 

The  case  —A  ——  2.405  corresponds  to  the  degenerate  case  of  the  first 
CO 


zero  of  the  zeroth  ordc 


Bessel  function 


Lorentz  force  and  averaged  radiation  power  gain  due 
to  the  first  harmonic.  Now,  we  can  deduce  that 


V\c  — — ^j-cri  AqAvq,  where  we  have  defined  Aq  = 
( K‘»  U°  ),  Vic  = ( Vlc  ),  Avd  = vD  - Vo  = 


-Wo 
Vd  - 


^ep 


V elc 


VO 


VeD  - VeO 


We  can  deduce  the  ac  gain  through  either  of  the 
two  equivalent  relations  (J\{ui)  — — V\c{u)/Cj\  = 

— real(f)i(a;))/(a>i/2).  Here  Vic(&)  and  V\{uj)  are  the 
first  harmonic  cosine  and  real  coefficient  of  the  Fourier 
expansion  of  the  limit  cycle  time  series  of  the  velocity. 
The  scaling  of  each  value  of  this  nondimensional  version 
nonlinear  ac  conductivity  to  a a whole  continuum  of  ex- 
perimental parameters  can  be  found  through  Tab.  1. 

Therefore,  given  that  there  are  resonances  in 
V\c(u),  cDi,  u>o)  vJjj  = c^o  and  that  the  peak  of  the 
I-V  curve  without  ac-drive  is  at  u)q  = 1 in  units  of  the 
geometric  mean  collision  frequency,  we  can  only  expect 
effects  of  the  first  resonance  on  the  I-V  characteristic 
in  the  form  of  additional  bumps  if  Q > 1.  Morover,  we 
qualitatively  deduce  that  if  V\c  is  slowly  varying  with  U)\ 
in  the  quasistatic  regime  O <<  1 that  one  can  expect  an 
approximate  linear  shift  in  vjj  with  increasing  E.g., 
a linear  shift  in  the  position  of  the  maximum  of  the  I-V 
curve  can  be  expected. 
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